Swaying threads of a solar filament 
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From recent high resolution observations obtained with the Swedish 1-m Solar 
Telescope in La Palma, we detect swaying motions of individual filament threads in 
the plane of the sky. The oscillatory character of these motions are comparable with os- 
cillatory Doppler signals obtained from corresponding filament threads. Simultaneous 
recordings of motions in the line-of-sight and in the plane of the sky give informa- 
tion about the orientation of the oscillatory plane. These oscillations are interpreted in 
the context of the magnetohydrodynamic theory. Kink magnetohydrodynamic waves 
supported by the thread body are proposed as an explanation of the observed thread 
oscillations. On the basis of this interpretation and by means of seismological argu- 
ments, we give an estimation of the thread Alfven speed and magnetic field strength 
by means of seismological arguments. 

Subject headings: Sun: filaments, Sun: oscillations 



1. Introduction 

When seen in a sub-arcs ec resolution, s olar prominences/filaments are resolved in to numerous 



threa d-like structures (e.g.. iLin et al]|2005|) . These threads are magnetic in nature (M artin et al 



2008|). Observations from the past 40 years reveal the ever present small-amplitude oscillations in 
prominences and filaments. T he observed oscillatory periods range from less than 1 minute up to a 
few hours (e.g., lEngvoldll2008l) . The fact that line-of-sight oscillations are observed in prominences 
beyond the limb and in filaments against the disk suggests that the planes of the oscillation may 
acquire various orientations relative to the local solar reference system. Recent high r esolution 



studie s show waves propagating along individual threads bo th in quiescent filaments (e.g.. LLin et al 



2007|) and in active region filaments (|Okamoto et al.ll2007|) . In the present study the authors aim 
to combine motions of thin threads along the line-of-sight, as derived from Doppler data, and 
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associated motions in the plane of the sky, which will yield information on the orientation of the 
plane of oscillation. 

Magnetohydrodynamic (MHD) waves have been extensively investigated and proposed as 
reasonable candidates t o explain detected period i cities in promin ences and filaments (the reader 
is referred to reviews by lOliver & B allesterl 12002c lBallesterll2006|) . By modeling a filament thread 
as a straight magnetic flux tube, partially filled wit h cool filament-li ke plas ma and embedded in 
a mu ch rarer and hotter coronal-like environment, iDfaz et al.l (|2002|) and iDymova & Ruderman 
(|2005|) studied the linear magnetoacoustic wave modes supported by such a structure in the (5 = 
appro ximation, where ft is the ratio of the gas pressure to the magnetic pressure. On the other 
hand. ISoler et al.1 (|2008|) considered a more simple model made of a homogeneous thread but took 
the {3 ± case into account and included additional effects such as nonadiabatic mechanisms and 
mass flow. Despite these differences in the modeling, all these works similarly conclude that, in 
the case of thin threads, the only trapped wave mode that is able to produce a significant nonax 



isymm etric, trans verse displacement of the flux tube is the so-called kink mode (|Edwin & Roberts 



1983b. Recently, Terrad as et al.l (|2008f) interpreted the oscillation of prominence threads reported 
by lOkamoto et al.l (|2007h as kink modes. The observations of filament thread oscillations from the 
Ha sequences reported in the present paper seem to be also consistent with a kink wave interpre- 
tation. Here, this possibility is explored and the implications of the observations in the context of 
the MHD theory are discussed. 

This paper is organized as follows. Section [2] describes the observations and the data reduc- 
tion. The observational results are presented in Section [3] Later, Section!?] contains a theoretical 
interpretation of the observed oscillations. Finally, the conclusions are given in Section |5] 



2. Observations and data reduction 



The target filament was observed with the Swedish 1-m Solar Telescope (SST. IScharmer et al 



2003|) on August 2, 2007, at solar coordinates S26W50 (see panel (a) of Figured)). At this helio- 
centric location of the target filament we observe the local horizontal and vertical directions under, 
respectively, 40° and 50° relative to line-of-sight. The observing field-of-view was centered on 
a relatively open section of the filament where threads were less densely packed and could be 
resolved individually (see panels (b) - (d) of FigureQ]). 

The Ha narrow-band filtergr ams were obtained from the Solar Optical Universal Polarimeter 
(SOUP, iTitle & Rosenberg! 1 198 lb tunable filter (FWHM 12.8 pm) of the Lockheed Martin Solar 
and Astrophysics Laboratory. After the SOUP prefilter (FWHM 8A) and before the SOUP filter, 
10% of the light was reflected by a beamsplitter to a phase-diversity camera pair. These cameras 
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Fig. 1. — (a): A sub-image extracted from a full-disk Ha image from the Catania Astrophysical 
Observatory. The white arrow indicates the position of the SST target. Panel (b) shows one SST 
image in Ha line center. Panels (c) and (d) are the corresponding sharpened Ha image and the 
Dopplergram. Ten swaying threads were selected from the Ha line center images. Three of them 
are marked by solid white lines in panels (c) and (d). The dashed white lines marked the locations 
of the cuts that were made across the swaying threads. 
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provided Ha broadband images which were used as an anchor channel in the post data reduction. 
The two phase-diversity cameras and the SOUP camera are Sarnoff lKxlK CCD cameras, which 
cover an 65x65 arc sec field of view. The three cameras were simultaneously exposed by means 
of an optical chopper. The fast Sarnoff cameras allow us to record 37 images per second at one 
line position. However, the time for tuning line positions in the SOUP is relatively long (about 8.4 
seconds). We therefore used two observing programs: (i) From 07:52UT to 08:30UT, the SOUP 
filter was alternating successively between three wavelength positions in the order of: Ha red wing 
(+0.3A), line core (6562. 8A), blue wing (-0.3A), line core, red wing and etc. The cadence for each 
wavelength is 18.5 seconds; (ii) From 08:32UT to 08:48UT, the SOUP filter was fixed in the Ha 
line center position, in order to observe filaments at high rate (cadence of 3.9 seconds). 



We applied the Multi-Object Multi-Frame Blind Deconvolution (MOMFBD: lvan Noort et al 



2005) image restoration technique to the two phase-diversity images and the SOUP images al- 
together for each SOUP cycle. Precise alignment of the sequentially recorded SOUP images is 
guaranteed through the phase-diversity channel in the MOMFBD restoration process, see e.g., 



Langangen et al.1 (120081) for details on the processing of a similar data set. This way, Ha Doppler- 
grams can be derived that are virtually free from misalignment errors due to seeing effects. The 
Dopplergrams were derived by fitting Ha line core and two wing intensities to a Gaussian profile. 



Beckers (1964) introduced the so-called cloud model to account for the influence of a chro- 



mospheric line on the line shift produced b y a separate ove rlying solar structure ("cloud"). A 



number of studies of this effect followed (cf. LMein et al.lll994|) . In a recent quantitative study Wiik 



and Engvold (2009, in preparation) simulate the impact of an artificial solar filament ("cloud") on 
an observed, high resolution chromospheric Ha line spectrum. The authors vary the parameters 
of the artificial filament, i.e. the optical thickness (r), line source function (S^), Doppler width of 
the line absorption (AA D ) and line-of-sight velocity (V/ os ). The study shows that unless the optical 
thickness t »1, the resulting line shift in the combined line profile is substantially less than what 
corresponds to the input line-of-sight velocity. In all cases of the investigated parameter ranges, for 
optical thicknesses in the range 0.5 - 1.0, the line shift derived from the resulting Ha line profiles 
yield generally too low line shift by about a factor 2. In the following we apply this correction 
factor to the derived line-of-sight velocities. 

We obtained sets of images in time series: Ha line core, red and blue wings with Ha broad- 
band images. By using the broadband images as reference, the restored Ha wing images and the 
line core images were aligned and de-stretched. The alignment between sequential images is to 
sub-pixel accuracy (i.e., better than 0.063 arcsec). We further applied a Laplacian kernel to the 
Ha images and the Dopplergrams with the aim to suppress the chromospheric background signals. 
Panels (b) - (d) of Figure \T\ show one Ha line center image and the corresponding sharpened Ha 
image and the Doppler image, respectively. 
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3. Observational results 

3.1. Swaying threads seen in Ha images 

The Ha filtergrams resolve many thin threads in the filament. These threads are more clearly 
seen in the sharpened Ha images. The Ha sequences with cadences of 3.9 s and 18.5 s show that 
some filament threads sway back and forth in the plane of sky. Ten such swaying threads were 
selected for investigation. Three of them are marked as solid white lines in panel (c) and (d) of 
Figure \T\ For each thread, two or three perpendicular cuts were made to measure the properties 
of possible waves (see panel (c) of same figure). The Ha time-slice diagram of each cut clearly 
shows the swaying signature, see panels (a)-(d) of Figure [21 

Filament threads undergo a variety of changes that make them hard to follow for more than 
5-10 minutes at a time. Their contrast c hanges and makes them appear and disappear within this 



time scale. They all move continuously (IZirker et all 1 1990 : Lin et al.ll2003|) presumably an effect 



of motion of their anchoring points in the photosphere. They also often mingle with neighboring 
threads and observers lose track of them. However, the swaying motion is a separate, characteristic 
feature of the threads, which may be monitored and studied within the limited time frame. 

The temporal variations of the positions of the threads are fitted by sine curves, see panel 
(c)-(f) of Figure |2l from which the period (P) and the amplitude (A) of the wave are derived. 
Assuming the wave is propagating along a such thread, the phase difference between the two fitted 
curves represents the time interval when the wave passes through the two cuts. We derive the time 
interval (T) from the maximum cross-correlation of the two curves. Taking one of the two curves 
as a reference, we shift the other curve according to the derived time interval. The deviations of 
corresponding points between the reference curve and the shifted curve give the uncertainty in the 
measurement of time interval. Given the distance (L) between the two cuts, one obtains the phase 
velocities of the waves from V ps = Jp. The uncertainties of the phase velocity can be calculated 
from the uncertainties in the time interval, since the distance (L) between two cuts is a constant. 



Tabled] lists the measured wave properties of the ten selected Ha swaying threads. These are 

ttA 
P 



short period waves (mean period ~ 3.6 min) and the mean velocities of the swaying motions (^fr) 



are about 2 km s 



-l 



3.2. Thread oscillations derived from Dopplergrams 

The 40-min Ha Dopplergrams with the cadence of 18.5 s allow us to look for the oscillations 
in the line-of-sight direction. Thread #1 - #3 are selected from this sequence. As two examples, 
Figure [3] shows the time-slice diagrams in both Ha and the Doppler for thread #2 and #3. 
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Fig. 2. — Cuts "c3" and "c4" are perpendicularly crossing thread #2. Cuts "c5" and "c6" are 
crossing thread #3, see panel (c) of Figured! The Ha time-slice diagrams of these four cuts are 
shown in panel (a) and (b). The swaying patterns of these two threads are marked by the dashed 
lines. As two examples, the swaying pattern seen in cuts "c4" and "c6" are further shown in 
panel (c) and (d) as the temporal variations of the positions of the two threads (plus signs). These 
variations are periodic and are fitted by sine curves (solid lines), from which the period and the 
amplitude of the wave can be derived. Panel (e) and (f) show the fitted sine curves of both cuts for 
each thread. The phase difference between the two corresponding cuts can be used to measure the 
phase velocities of the waves, see text for details. 
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Fig. 3. — Panels (a) and (c): The time-slice diagrams of thread #2 and #3 in Ha line center, 
sharpened Ha line center images and the Dopplergrams. The two short sections of time- slice 
diagrams in each panel correspond to the period where the swaying effects of the threads are visible 
in Ha line center (c.f., Figure [2]). Due to the swaying effect, these two sections were reconstructed 
by adjusting the thread positions slightly. Panels (b) and (d): The time variation of the mean 
Doppler velocities averaged over each slice and the corresponding power spectrum. The dashed 
lines show the 90% confidence level. The detected periods are 5.3 minutes for thread #2 and 4.7 
minutes for thread #3. 
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Because the position of a slice (marked as some white solid lines in panel (c) of Figured)) is 
defined from one image and is fixed, when plotting a time-slice diagram, the filament thread may 
move out of such a fixed slice "window" due to the swaying effect. Therefore, we reconstruct the 
short sections of the time-slice diagrams, where the swaying threads are seen in Ha, by adjusting 
the slice positions slightly (see panel (a) and (c) of Figure [3]). The corresponding short section of 
the Doppler time-slice diagrams are reconstructed accordingly. 

For thread #2 and #3, the temporal variations of their mean Doppler velocities averaged over 
each slice indicate that the threads oscillate in the line-of-sight (see panels (b) and (d) of Figured]). 
From the power spectrum, the waves are short period (~ 5 min). Compared with thread #2 and #3, 
the Doppler signals of thread #1, which is not shown here, are relatively weak. 

The amplitude of the derived line-of-sight velocity ranges typically between 1 and 2 kms" 1 . 
For threads #2 and #3, they are 0.8 kms" 1 and 1.9 kms -1 , respectively. The quantitative study of 
Wiik and Engvold (2009, in preparation) also demonstrates that in cases when the optical thick- 
ness of the artificial filament r ~ 1 the velocity fluctuation caused by the fine structure of the 
chromosphere below is suppressed to yield an uncertainty cr x 0.5 kms" 1 . Applying this intrinsic 
error (0.5 kms" 1 ) caused by the chromosphere below and the "cloud model"-correction, the two 
line-of-sight velocities raise to 1.6+0.5 kms" 1 and 3.8+0.5 kms" 1 . 



3.3. Combination of the observed oscillations in the two orthogonal directions 

The two selected Ha threads are found oscillating both in the plane of sky and in the line- 
of-sight, with a similar period. Since the location of the target filament is S26W50, the measured 
swaying motion and oscillation of velocity in line-of-sight represent inevitably the projected com- 
ponents of the oscillating motions. As noted earlier, the oscillations of the filament threads are 
probably polarized and the plane of oscillation may attain various orientation relative to the local 
reference system. Swaying motions will be most clearly observed when a thread sways in the plane 
of the sky while Doppler signals will be strongest for oscillations along the line of sight. 

Combining the corresponding components derived from the swaying motion of the two threads, 
which are 1.8+0.2 kms" 1 and 2.3+0.2 kms" 1 , and the line-of-sight oscillatory motion, 1.6+0.5 
kms" 1 and 3.8+0.5 kms" 1 , we may derive the full vectors as 2.4+0.5 kms" 1 and 4.4+0.5 kms" 1 
which are directed at angles 42+10° and 59+4°, respectively, relative to the plane of the sky. Since 
the heliocentric position of the target filament is W50, the corresponding local vertical direction 
becomes 50° relative to line-of-sight and we conclude that these two threads oscillate in planes 
which are reasonably close to vertical. It is, however, not possible based on only two cases to draw 
any general conclusion about the orientation of planes of oscillation of filament threads. 
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From Table Q], one notices that the amplitudes of the swaying motions are very small. This 
calls for very high spatial resolution in observations, which may explain why it has not been no- 
ticed in earlier studies. On the other hand, for some threads, the amplitudes of the waves passing 
through two cuts are notably different. In other words, the amplitudes of the waves are changing 
while they propagate along these threads. Apparent changes can be due to spatial damping of 
the waves in addition to the noise in the data. The damping phenomenon appears very common 
in solar filaments/prominences with a damping t ime corresponding to from 1 to 4 periods (e.g., 



Molownv-Horas et al. 



199c 



Terradas et al. 



Terradas et al. 



l200ll 



2002) . Various possible theoretical mechanis ms have 



2002: ISoler et al.ll2008l : lArregui et all l2008bl) . among 



been proposed (e.g., 

which the mechanis m of resonant absorp t ion seems to be th e most efficient one for the damping 
of kink waves (e.g., lArregui et al.l 120081* ISoler et al.l 120091) . It is also expected this mechanism 
to be the dominant one for the spatial damping of propagating kink waves, although a theoretical 
investigation on this issue is needed. Alternatively, the observed change of the wave amplitude 
might also be caused by a rotation of the oscillating planes. We will look for such evidence from 
higher temporal resolution data (cadence less than 2 sec) in a following filament seismology study. 



4. Theoretical interpretation of the oscillations 
4.1. Kink magnetohydrodynamic waves in a magnetic flux tube 



Next, in order to perform an interpretation of the observed oscillations in terms of linear 
magnetohydrodynamic waves, we consider a simple model representing a prominence thread. We 
assume a configuration made of a straight and homogeneous cylinder of length L and radius a, 
filled with filament-like plasma and embedded in an also homogeneous coronal environment. The 
magnetic field is taken homogeneous and orientated along the tube axis (see a scheme of the 
model in Fig H]). In this model, we neglect the longitudinal inhomogeneity of the thread plasma 
for simplicity and consider /3 = 0, so magnetic effects are dominant over gas pressure effects. By 
using cylindrical coordinates, we assume that the density profile, p , only depends on the radial 
coordinate, r, as follows, 

Pi, if r < a, 
p e , if r > a, 

where p, and p e are the internal and external densities respectively. 



Po if) = 



(1) 



The linear magnetoh ydrodynamic wave modes supported by this flux tube model have been 
investigated in detail (e.g.. ISpruitll 19821 : Edwin & Robertsll 19831 : ICallyll 1986|) . Among the possible 
wave modes, the kink mode is the only one that can produce a significant transverse displacement 
of the cylinder axis, and so it is the best candidate for an interpretation of the present oscillations. 
Moreover, the kink mode produces short-period oscillations of the order of minutes for typical fil- 
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Table 1: For the ten selected Ha swaying threads, the periods (P), phase velocities (V p/! ) and the 
amplitudes (A) of the possible waves were measured. Thread #l-#3, indicated in Figure [3] (c), are 
selected from the time series with the cadence of 18.5 seconds. Thread #4-#10 are from the time 
series with the cadence of 3.9 seconds. 
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Fig. 4. — Sketch of the model configuration. 



ament plasma conditions, which is also compatible with the measured periods. Since the observed 
threads appear to be very thin and long structures, it seems reasonable to consider the so-called 
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Fig. 5. — a) Ratio c\lv 2 Ai (solid line) as a function of the density contrast, c. The dotted line 
corresponds to the value of the ratio c 2 k lv 2 Ai for c — » oo. b) Magnetic field strength, B , as a function 
of the internal density, p h corresponding to some selected threads: thread 1 (solid line), thread 3 
(dotted line), thread 5 (dashed line), and thread 7 (dash-dotted line). 

thin tube approximation, i.e., L » a. In such an approximation, the kink mode phase velocity 
(also called the kink speed), namely q, is given by 



Ck = 



2BI 



n(Pi+p e y 



(2) 



where B Q is the magnetic field strength and /u = An x 10~ 7 A^A" 2 is the magnetic permittivity. 
Equation © can be rewritten in terms of the internal Alfven speed, va; = Bq/ y/jJpi, and the density 
contrast, c = p,/p e , as follows 



Ck = v A ; 



2c 

cTT' 



(3) 



Figured displays the ratio c 2 k lv 2 Ai as a function of c. We see that for large density contrast the 
curve becomes flat and so the ratio c 2 k lv 2 Ai is then almost independent of the density contrast. For 

typical coronal and filament densities one has c - 200. In such a situation, the factor ^ * V2 
in Equation © and so the kink speed directly depends on the internal Alfven speed, 



Ck 



(4) 



Now, assuming that thread oscillations observed from the Ha sequences are the result of a 
propagating kink mode, the measured phase velocities, V p h, can be related to the kink speed, so 
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Ck = Vph- Therefore, it is possible to use Equation © to give a seismological estimation of the 
thread Alfven speed, 



VAi 



V ph /^2. 



(5) 



The use of the MHD seismology to obtain information of the plasma physica l conditions has been 
applied by some authors in the context of coronal lo op oscillations (e.g., iNakariakov & Ofman 



2001 



Arregui et al. Il200l 2008al : lGoossens et al.ll2008b and prominence oscillations (|Terradas et all 



2008J). Since the density contrast in the case of coronal loops typically has a small value, the factor 
with the density contrast cannot be dropped from equation © an d so it is not possible to give 
a direct estimation of the Alfven speed from the observed kink speed. However, the much larger 
density of filament threads allows us to provide a more accurate determination of the Alfven speed. 



4.2. Determination of the Alfven speed and the magnetic field strength 

Table [2] contains an estimation of the Alfven speed of each swaying thread observed in the 
Ha sequences by taking expression © into account. We see that the Alfven speed varies in a 
wide range, which suggests that physical properties (i.e., density, magnetic field strength, etc.) 
significantly changes in different threads of the filament. Once the Alfven speed is determined, 
we have a relation between the magnetic field strength and the thread density, i.e., B = v A , yjjjpi- 
Since the thread density is an unknown parameter, we cannot uniquely determine the magnetic 
field strength or vice versa. Figure [5h> shows the dependence B with p, corresponding to four 
selected threads. The density has been varied in a wide range, 10~ 12 kg mr 3 < p t < 10~ 9 kg m -3 , 
and subsequently the magnetic field strength is in the range 0.1 G < B < 20 G. This range agrees 
with the usually reported values of the magnetic field strength in prominences and filaments (e.g., 



Patsourakos & Vial 2002|) . From Figure it can be also observed that the larger the density, the 



larger the magnetic field strength, which is compatible with the idea that active region filaments 
have stronger magnetic fields. By assuming a typical value of p, ■ = 5 x 10" 11 kg m~ 3 for the thread 
density and considering the results for the ten studied threads, one obtains a magnetic field strength 
which varies in the range 0.5 G < B < 5 G, see Tabled Although we have to be cautious with the 
values obtained in this very rough estimation, the present results suggest that the magnetic field in 
the filament is not homogeneous and could largely vary between threads. Also, this could happen 
for the density, but the results are more sensible to variations of the magnetic field strength than 
variations of the density. 
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va,- [kms" 1 ] 
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1.1+0.2 


5 


41+6 


3.2+0.7 


10 


20+6 


1.6+0.7 



Table 2: Estimation of the Alfven speed, va,, and magnetic field strength, B , of threads listed 
in Table I by considering expression ©. The value of B has been obtained by assuming p, = 
5x 10" 11 kgm" 3 . 



5. Concluding remarks 

In this study, the SST Ha intensity images and the Dopplergrams reveal many thin threads of 
a quiescent filament. The Ha line center images in time series show a number of threads swaying 
back and forth in the plane of sky, with a mean period of 3.6 minutes, a mean phase velocity of 33 
km s~ ! and a mean amplitude velocity of 2 kms" 1 . A si milar swaying p a ttern o f fila ment threads 



is also seen in a recent Hinode Ca II H observation by lOkamoto et al.l (|2007|) and iBerger et al. 



(2008). Okamoto et al. found some horizontal threads in an active region prominence oscillating 



vertically in the dark plane of sky, with periods of 130-250 s and amplitudes from 400 to 1800 km. 
Their measured periods are consistent with the current study. However, thei r amplitudes appear 



larger which might be because their prominence target is in an active region. IBerger et alj (|2008) 
also found transverse oscillations in quiescent prominences with intermediate periods of 20-40 
minutes, propagating speed of 10 km s _1 and amplitudes of 100-250 km. The longer duration of 
their observations (several hours) favors the detection of longer periods. It should be no t ed tha t 



similar swaying motions have also been observed in solar spicules by iDe Pontieu et al] (|2007l) . 
who interpreted the transverse motions as Alfvenic waves. 

Two filament threads in the current dataset are found both swaying in the plane of sky and 
oscillating in the line of sight with similar periods. Assuming they are the two projected compo- 
nents of the oscillating motions, the orientation of the oscillating plane can be derived. This will 
be further investigated in a follow-up study based on a higher temporal resolution data (cadence 
less than 2 s). 

We have interpreted the observed thread oscillations in the Ha sequences in terms of kink 
magnetohydrodynamic wave modes propagating along the threads. In the context of this interpre- 
tation, we have provided a determination of the Alfven speed of the studied threads. Subsequently, 
we have assumed a typical value for the plasma density and have estimated the magnetic field 
strength in the filament. High resolution observations of filament threads suggest that the plasma 
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parameters may be varying along individual threads, which may affect, somehow, the derived mag- 
netic field strength. 
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